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We calculate the fully differential next-to-next-to-leading-order (NNLO) corrections to vector- 
boson fusion (VBF) Higgs production at proton colliders, in the limit in which there is no cross¬ 
talk between the hadronic systems associated with the two protons. We achieve this using a new 
“projection-to-Born” method that combines an inclusive NNLO calculation in the structure-function 
approach and a suitably factorised next-to-leading-order (NLO) VBF Higgs plus 3-jet calculation, 
using appropriate Higgs plus 2-parton counterevents. An earlier calculation of the fully inclusive 
cross section had found small NNLO corrections, at the 1% level. In contrast, the cross section 
after typical experimental VBF cuts receives NNLO contributions of about 5-6%, while differential 
distributions show corrections of up to 10-12% for some standard observables. The corrections are 
often outside the NLO scale-uncertainty band. 


Following the discovery in 2012 of the Higgs boson [1] 
[2] , one of the main tasks for particle physics today is the 
accurate determination of its properties and couplings. 
For the coming decade at least, these studies will take 
place at CERN’s Large Hadron Gollider (LHC). 

The most relevant production channels for the Higgs 
boson at the LHC are gluon fusion (ggH), vector-boson 
fusion (VBFH), production in association with a vector 
boson (VH) and with a top-quark pair (ttH) [5]. VBFH 
is special for a number of reasons HIS]: it has the largest 
cross section of the processes that involves tree-level pro¬ 
duction of the Higgs boson (and is second largest among 
all processes); it has a distinctive signature of two for¬ 
ward jets, which makes it possible to tag the events and 
so identify Higgs decays that normally have large back¬ 
grounds, e.g. H —>■ T+r“; the Higgs transverse momen¬ 
tum is non-zero even at lowest order, which facilitates 
searches for invisible decay modes PEI; and it also brings 
particular sensitivity to the charge-parity properties of 
the Higgs boson, and non-standard Higgs interactions, 
through the angular correlations of the forward jets [S] . 

Given the key role of VBF Higgs-boson production at 
the LHC, it is of paramount importance to have a pre¬ 
cise prediction for its production. The total cross sec¬ 
tion was calculated to NNLO in Refs. mm using the 
structure-function approach m, showing small correc¬ 
tions relative to the NLO and tiny renormalisation and 
factorisation scale uncertainties, well below 1%. However 
experimental measurements are necessarily restricted to 
a subset of phase space. In particular, because of their 
use of transverse-momentum cuts on the forward tagging 
jets, one might imagine that there are important NNLO 
corrections, associated with those jet cuts, that would 
not be seen in a fully inclusive calculation. Currently, 
the fully differential VBFH cross section is known only 
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to NLO [12]. It appears to have small scale uncertainties. 

One can think of the VBFH process, represented at the 
Born level in Fig. as involving two Deeply Inelastic 
Scatterings (DIS), one for each of the incoming protons. 
Each DIS process produces a vector boson, or Z, 
and the fusion of the vector bosons produces a Higgs bo¬ 
son. The structure-function approach m used for the 
NNLO total cross section liiini assumes that the up¬ 
per and lower hadronic sectors factorise from each other, 
i.e. that there is no cross-talk between them. Factorisa¬ 
tion is believed to be accurate to better than 1% in the 
experimentally relevant kinematic region 

The reason that the structure function approach does 
not provide a fully differential cross section is related to 
the fact that the DIS coefficient functions used in the 
calculation implicitly integrate over hadronic final states. 
In this letter, we introduce a new “projection-to-Born” 
approach to eliminate this limitation, and thus provide 
the first fully differential NNLO calculation of VBF Higgs 
production in the factorised approximation. 

Let us start by recalling that the cross section in the 
structure-function approach is expressed [SHII] as a sum 
of terms involving products of structure functions, e.g. 
F 2 {xi,Ql)F 2 {x 2 ,Ql), where Qf = -qf > 0 is given in 
terms of the 4-momentum Qi of the (outgoing) exchanged 
vector boson i. The Xi values are fixed by the relation 
Xi = —Q‘ll{2Pi.qi)^ where Pi is the momentum of proton 
i. To obtain the total cross section, one integrates over all 
Qi, q 2 that can lead to the production of a Higgs boson. 
If the underlying upper (lower) scattering is Born-like, 


^ The factorisation of the two sectors is exact if one imagines two 
copies of QCD, QCDi and QCD 2 , respectively for the upper 
and lower sectors, where each of the two QCD copies interacts 
with the electroweak (EW) sector, but not with the other QCD 
copy. This observation could be exploited, for example, in auto¬ 
mated calculations and for determining corrections beyond the 
factorised approximation. 
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(b) NNLO "inclusive" part (from structure function method) 


two loop 


projeeted double real 


(c) NNLO "exclusive" part (from VBF H+3j@NLO) 


double real double-real counterevent 



FIG. 1: (a) Illustration of the Born VBFH process, (b) NNLO corrections to the upper sector of the VBF process, from the 
“inclusive” part of our calculation, (c) Corresponding “exclusive” part. The double-real and one-loop single-real counterevents 
in the exclusive part cancel the projected double-real and one-loop single-real contributions in the inclusive part. In the 
“projected” and “counterevent” contributions, the dashed lines corresponds to the full set of parton momenta that are integrated 
over (for the structure functions, this integral is implicit in the derivation of the coefficient functions), while the solid lines 
correspond to the partons that are left over after projection to Born-like kinematics and then passed to the analysis. The 
projection does not change the direction of initial partons and so the corresponding incoming dashed lines are implicit. 


quark —>■ quark -1- V, then it is straightforward to show 
that knowledge of the vector-boson momentum qi (^ 2 ) 
uniquely determines the momenta of both the incoming 
and outgoing (on-shell) quarks, 

Pin,i — Pout,j — Qi • (1) 

We exploit this feature in order to assemble a full cal¬ 
culation from two separate ingredients. For the first one, 
the “inclusive” ingredient, we remain within the struc¬ 
ture function approach, and for each set of qi and <72 use 
Eq. Q to assign VBF Born-like kinematics to the up¬ 
per and lower sectors. This is represented in Fig. & 
(showing, for brevity, just the upper sector): for the 
two-loop contribution, the Born kinematics that we as¬ 
sign corresponds to that of the actual diagrams; for the 
tree-level double-real and one-loop single-real diagrams, 
it corresponds to a projection from the true kinematics 
(2 —>■ iL -I- n for n = 3,4) down to the Born kinemat¬ 
ics (2 —>• H + 2). The projected momenta are used to 
obtain the “inclusive” contribution to differential cross 
sections. Note that the Higgs momentum is unaffected 
by the projection. 

Our second, “exclusive”, ingredient starts from the 
NLO fully differential calculation of vector-boson fusion 
Higgs production with three jets miiii], as obtained in 
a factorised approximation, i.e. where there is no cross¬ 
talk between upper and lower sectors.^ Thus each par- 
ton can be uniquely assigned to one of the upper or lower 
sectors and the two vector-boson momenta can be unam¬ 
biguously determined. For each event in a Monte Carlo 
integration over phase space, with weight w, we add a 


^ The NLO calculation without this approximation is given in 
Ref. [18]. 


counterevent, with weight —w, to which we assign pro¬ 
jected Born VBF kinematics based on the vector-boson 
momenta and Eq. Q. This is illustrated in Fig. [^. 
From the original events, we thus obtain the full mo¬ 
mentum structure for tree-level double-real and one-loop 
single-real contributions. Meanwhile, after integration 
over phase space, the counterevents exactly cancel the 
projected tree-level double-real and one-loop single-real 
contributions from the inclusive part of the calculation. 
Thus the sum of the inclusive and exclusive parts gives 
the complete differential NNLO VBFH result.^ 

For the implementation of the inclusive part of the cal¬ 
culation, we have taken the phase space from POWHEG’s 
Higgs plus two-jet VBF calculation [20], while the matrix 
element has been coded with structure functions evalu¬ 
ated using parametrised versions im |22] of the NNLO 
DIS coefficient functions [231 - 125] integrated with MOPPET 
vl.1.5 [55]. We have tested our implementation against 
the results of one of the codes used in Ref. imnj and 
found agreement, both for the structure functions and the 
final cross sections. We have also checked that switching 
to the exact DIS coefficient functions has a negligible im¬ 
pact. A further successful comparison of the evaluation of 
structure functions was made against APPEL v.2.4. I|27|. 

For the exclusive part of the calculation, as a starting 
point we took the NLO (i.e. fixed-order, but not parton- 
shower) part of the POWHEG iL-|-3-jet VBF code [17], it¬ 
self based on the calculation of Ref. m, with tree-level 
matrix elements from MadGraph 4 [55]. This code al¬ 
ready uses a factorised approximation for the matrix ele¬ 
ment, however for a given phase-space point it sums over 


® Our approach can be contrasted with the differential NNLO 
structure-function type calculation for single-top production m 
in that we do not need any fully differential ingredients at NNLO. 
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matrix-element weights for the assignments of partons 
to upper and lower sectors. We therefore re-engineered 
the code so that for each set of 4-momenta, weights are 
decomposed into the contributions for each of the dif¬ 
ferent possible sets of assignments of partons to the two 
sectors. For every element of this decomposition it is 
then possible to unambiguously obtain the vector-boson 
momenta and so correctly generate a counterevent. The 
PQWHEG-BOX’s [m EO] “tagging” facility was particularly 
useful in this respect, notably for the NLO subtraction 
terms. To check the correctness of the assignment to 
sectors, we verified that as the rapidity separation be¬ 
tween the two leading jets increases, there was a decreas¬ 
ing relative fraction of the cross section for which partons 
assigned to the upper (lower) sector were found in the ra¬ 
pidity region associated with the lower (upper) leading 
jet. We also tested that the sum of inclusive and exclu¬ 
sive contributions at NLO agrees with the PQWHEG NLO 
implementation of the VBF H+2-iet process. 

To investigate the phenomenological consequences of 
the NNLO corrections, we study 13 TeV proton-proton 
collisions. We use a diagonal CKM matrix, full Breit- 
Wigners for the W, Z and the narrow-width approxima¬ 
tion for the Higgs boson. We take NNPDF 3.0 parton 
distribution functions at NNLO with as{Mz) = 0.118 
(NNPDF30_nnlo_as_0118) [ST], also for our LO and NLO 
results. We have five light flavours and ignore contribu¬ 
tions with top-quarks in the final state or internal lines. 
We set the Higgs mass to Mh = 125 GeV, compati¬ 
ble with the experimentally measured value |32]- Elec- 
troweak parameters are set according to known exper¬ 
imental values and tree-level electroweak relations. As 
inputs we use Mw = 80.398 GeV, Mz = 91.1876 GeV 
and Gp = 1.16637 x 10“® GeV“^. For the widths 
of the vector bosons we use F^ = 2.141 GeV and 
Tz = 2.4952 GeV. 

Some care is needed with the renormalisation and fac¬ 
torisation scale choice. A natural option would be to use 
Qi and Q 2 as our central values for the upper and lower 
sectors, respectively. While this is straightforward in the 
inclusive code, in the exclusive code we had the limitation 
that the underlying PDWHEG-BOX code can presently only 
easily assign a single scale (or set of scales) to a given 
event. However, for each PQWHEG phase-space point, we 
have multiple upper/lower classifications of the partons, 
leading to several {Qi, < 52 } pairs for each event. Thus the 
use of Qi and Q 2 would require some further degree of 
modification of the POWHEG-BDX, which we leave to future 
work. We instead choose a central scale that depends on 
the Higgs transverse momentum pt,H- 



This choice of po is usually close to 'jQxQi- It represents 
a good compromise between satisfying the requirement of 
a single scale for each event, while dynamically adapting 
to the structure of the event. In order to estimate missing 


^(no cuts) ^(VBF cuts) jp^j 


LO 

4 0*^2 

^.UOZ _o.069 

0.957 

NLO 

_Q 023 

0.876 

NNLO 

Q QQO +0.016 
0.000 -0.012 

0.826 


TABLE I: Cross sections at LO, NLO and NNLO for VBF 
Higgs production, fully inclusively and with VBF cuts. The 
quoted uncertainties correspond to scale dependence, while 
statistical errors at NNLO are about 0.1% with VBF cuts 
and much smaller without. 


higher-order uncertainties, we vary the renormalisation 
and factorisation scales symmetrically (i.e. keeping pp = 
Pp) by a factor 2 up and down around 

To pass our VBF selection cuts, events should have at 
least two jets with transverse momentum pt > 25 GeV; 
the two hardest (i.e. highest pt ) jets should have absolute 
rapidity |?/| < 4.5, be separated by a rapidity > 

4.5, have a dijet invariant mass > 600 GeV and 

be in opposite hemispheres < 0). Jets are de¬ 

fined using the anti-fct algorithm |33) . as implemented in 
Fast Jet v3.1.2 [31], with radius parameter R = 0.4. 

Results are shown in table [T| for the fully inclusive cross 
section and with our VBF cuts. One sees that the NNLO 
corrections modify the fully inclusive cross section only 
at the percent level, which is compatible with the find¬ 
ings of Ref. [5|. However, after VBF cuts, the NNLO 
corrections are about 5 times larger, reducing the cross 
section by 5 — 6% relative to NLO. The magnitude of the 
NNLO effects after cuts implies that it will be essential 
to take them into account for future precision studies. 
Note that in both the inclusive and VBF-cut cases, the 
NNLO contributions are larger than would be expected 
from NLO scale variation. 

Differential cross sections are shown in Fig. for 
events that pass the VBF cuts. From left to right, the 
plot shows the transverse momentum distributions for 
the two leading jets, ptj^ and Pt,j 2 ^ Higgs boson, 

Pt,H, and the distribution for the rapidity separation be¬ 
tween the two leading jets, Ayj^ j^. The bands and the 
patterned boxes denote the scale uncertainties, while the 
vertical error-bars denote the statistical uncertainty. The 
effect of the NNLO corrections on the jets appears to be 
to reduce their transverse momentum, leading to nega¬ 
tive (positive) corrections in regions of falling (rising) jet 
spectra. One can see effects of up to 10 — 12%. Turn¬ 
ing to pt,H, one might initially be surprised that such an 
inclusive observable should also have substantial NNLO 
corrections, of about 8% for low and moderate pt.H- Our 


^ We verified that an expanded scale variation, allowing pp ^ pp 
with I < pp/pp < 2, led only to very small changes in the 
NNLO scale uncertainties for the VBF-cut cross section and the 
Pt,H distribution. 








4 


do/dptj^ [pb/GeV] 


do/dpt [pb/GeV] 


do/dpt,H [pb/GeV] 


10-2 

10 -^ 

10-4 

1.1 


0.9 

0.8 

50 100 150 200 250 300 20 40 60 80 100 120 140 160 0 50 100 150 200 250 300 

Pt.j, [GeV] Pt [GeV] [GeV] 




. .. 

4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 




FIG. 2: From left to right, differential cross sections for the transverse momentum distributions for the two leading jets, Pt,ji 
and Pt,j 2 i for the Higgs boson, pt,H, and the distribution for the rapidity separation between the two leading jets, 


interpretation is that since NNLO effects redistribute jets 
from higher to lower pt’s (cf. the plots for pt,ji and Pt,j^), 
they reduce the cross section for any observable defined 
with VBF cuts. As pt^H grows larger, the forward jets 
tend naturally to get harder and so automatically pass 
the pt thresholds, reducing the impact of NNLO terms. 

As observed above for the total cross section with VBF 
cuts, the NNLO differential corrections are sizeable and 
often outside the uncertainty band suggested by NLO 
scale variation. One reason for this might be that NLO 
is the first order where the non-inclusiveness of the jet 
definition matters, e.g. radiation outside the cone modi¬ 
fies the cross section. Thus NLO is, in effect, a leading- 
order calculation for the exclusive corrections, with all 
associated limitations. 

To further understand the size of the NNLO correc¬ 
tions, it is instructive to examine a NLO plus parton 
shower (NLOPS) calculation, since the parton shower 
will include some approximation of the NNLO correc¬ 
tions. For this purpose we have used the POWHEG VBF 
iL-|-2-jet calculation [7D], showered with PYTHIA version 
6.428 with the Perugia 2012 tune [3S]. The POWHEG part 
of this NLOPS calculation uses the same PDF, scale 
choices and electroweak parameters as our full NNLO 
calculation. The NLOPS results are included in Fig. 
at parton level, with multi-parton interactions (MPI) 
switched off. They differ from the NLO by an amount 
that is of a similar order of magnitude to the NNLO 
effects. This lends support to our interpretation that fi¬ 
nal (and initial)-state radiation from the hard partons 
is responsible for a substantial part of the NNLO correc¬ 
tions. However, while the NLOPS calculation reproduces 
the shape of the NNLO corrections for some observables 


(especially pt,H), there are others for which this is not 
the case, the most striking being perhaps Par- 

ton shower effects were also studied in Ref. [35], using 
the MC@NLO approach EZj. Various parton showers 
differed there by up to about 10%. 

In addition to the NNLO contributions, precise phe¬ 
nomenological studies require the inclusion of EW con¬ 
tributions and non-perturbative hadronisation and MPI 
corrections. The former are of the same order of magni¬ 
tude as our NNLO corrections [13]. Using Pythia 6.428 
and Pythia 8.185 we find that hadronisation corrections 
are between —2 and 0%, while MPI brings up to +5% 
at low Pt’s. The small hadronisation corrections appear 
to be due to a partial cancellation between shifts in pt 
and rapidity. We leave a combined study of all effects 
to future work. The code for our calculation will also be 
made public. 

With the calculation presented in this letter, differ¬ 
ential VBF Higgs production has been brought to the 
same NNLO level of accuracy that has been available for 
some time now for the ggH [SSJ [33] and VH [3D] pro¬ 
duction channels. This constitutes the first fully differ¬ 
ential NNLO 2 —3 hadron-collider calculation, an ad¬ 
vance made possible thanks to the factorisable nature of 
the process. The NNLO corrections are non-negligible, 
5-10%, i.e. an order of magnitude larger than the cor¬ 
rections to the inclusive cross section. Their size might 
even motivate a calculation one order higher, to N^LO, 
to match the precision achieved recently for the ggH to¬ 
tal cross section [41]. With the new “projection-to-Born” 
approach introduced here, we believe that this is within 
reach. It would also be of interest to obtain NNLO plus 
parton shower predictions, again matching the accuracy 
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achieved recently in ggH |321133] . 
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